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Abstract. We present results from a series of observational tests to 3D and ID solar models. 
In particular, emphasis is given to the line formation of atomic oxygen lines, used to derive 
the much debated solar oxygen photospheric abundance. Using high-quality observations 
obtained with the Swedish Solar Telescope (SST) we study the centre-to-limb variation of 
the O i lines, testing the models and line formation (LTE and non-LTE). For the O i 777 nm 
triplet, the centre-to-limb variation sets strong constraints in the non-LTE line formation, 
and is used to derive an empirical correction factor (Sh) to the classical Drawin recipe for 
neutral hydrogen collisions. Taking advantage of the spatially-resolved character of the SST 
data, an additional framework for testing the 3D model and line formation is also studied. 
From the tests we confirm that the employed 3D model is realistic and its predictions agree 
very well with the observations. 
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1. Introduction 

Three-dimensional, time dependent hydrody- 
namical simulations of stellar atmospheres rep- 
resent a paradigm c hange in the m odeling of 
stellar atmospheres dAsplunol I2005I) . Unlike 
their one-dimensional counterparts, they treat 
convection self-consistently without the need 
for broadening parameters such as micro and 
macroturbulence. Their realistic description of 
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the motions and velocities in the photosphere 
leads to an excellent agreement between the 
predicted and observ ed line shapes and shifts 
(|AsplundetaTll2000h. 

lAsplund etai] (l2005h have revised the so- 
lar chemical composition, using photospheric 
lines t he 3D solar model of IStein & Nordlundl 
(1998). This revision implied a lower metal 
content, in particular C, N and O, which 
caused significant in the agreement of solar 
interior models with results f rom helioseis- 
mology dBahcall et all 120051: iBasu & Antial 
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20081). The n e w ch emical composition of 



Asplu nd et aD (120091) results in a slightly 



higher solar metalicity, but per se is not enough 
to reconcile the solar interior models with he- 
lioseismology. To trust the results from the 3D 
models, more testing is needed. Our objective 
is to te st the 3D model used by Asplu nd et al.1 
(2009), in particular for the oxygen lines, de- 
terminant in the downward revision of the 
solar abundances. This contribution su mma- 
rizes the findings of iPereira et alJ (l2009bl) and 
iPereiraet all J2009a). where the centre-to-limb 
variation and spatially-resolved line formation 
of oxygen lines is tested. 

2. Observations and photospheric 
models 

Using the Swedish 1-m Solar Telescope (SST) 
with the TRIPPEL spectrograph, we obtained 
high spatial and spectral resolution observa- 
tions at several positions in the solar disk. The 
lines observed include five atomic oxygen lines 
at 615.81 nm, 630.03 nm and the three lines 
around 777 nm. Stray light in the telescope was 
estimated and corrected for . The o bservations 
are detailed in lPereira et al.l (l2009bl) . 

We test these observations against a new 
3D model atmosphere (Trampedach et al. 
2009, in preparation). I t was computed with the 
stagger MHD code dNordlund & Gals g aard , 
119951: iGudiksen & Nordlundl l2005h . and con- 
sists of a 240 3 grid with a physical size of 
6x6x4 Mm. It includes an improved treatment 
of radiation, with a 12-bin multi-group opacity 
binning scheme. For the line formation calcu- 
lations the original simulation was interpolated 
to a 50x50x82 grid to save computing time. 
The simulation snapshots used here cover «45 
minutes of solar time. 

In addition, we also compare the centre-to- 
limb variation predicted b y two ID models: the 
semi-e mpirical model of Holweger & Miiller 
(11974 and the LTE, line-bl anketed marcs 
model dGustafsson et all l2008h . 

Oxygen line formation calculations were 
done in LTE and NLTE, u sing our LTE 
code and the mult i3d 
Botnen & Carlssorl 



code jBotnenl 11997 
Il999t lAsplund etal 



20031) . A 23-level model atom was used (see 



details in IPereira et all l2009al) . For the Oi 
615.81 nm and [Oi] 630.03 nm lines we show 
the LTE results, because the NLTE effects are 
weak in those lines. For the Oi 777 nm lines 
NLTE calculations were carried out for several 
values of S h, a scaling factor for the col lisions 
with Hi from the classical fo r mulae (iDrawini 
119681: ISteenbock & Holwegerl Il984l) . Nearby 
blends were included in the calculations of 
the 615.81 nm line (in particular, CN and C2 
molecular blends) and the 630.03 nm line (in 
particular, the Ni 1 blend). 

3. Centre-to-limb variation 

In Fig. Q] we show the centre-to-limb varia- 
tion of some of the lines. For the 615.81 nm 
and 630.03 nm lines the 3D model predictions 
agree well with the observations. For these 
lines there is only a small difference between 
the different models. 

For the Oi 777 nm lines, the efficiency 
of the collisions with neutral hydrogen in 
the NLTE calculations is u nknown. Following 
lAllende Prieto et al~ (2004), we use the centre- 
to-limb variation of these lines to derive an em- 
pirical estimate of Sh, the multiplier factor by 
the classical collision rates. We show only one 
representative line, and only for the 3D model. 
Of the Sh values used, we find that Sh — 1 
gi ves the best agreemen t with the observations. 
In lPereira et al.1 (l2009al) a further refinement of 
this value is made, and the best fitting value is 
found to be S h ~ 0.85. For these lines, one 
can very clearly rule out LTE as an accept- 
able approximation. The 777 nm results for the 
ID Holweger-Miiller model are very similar to 
those of the 3D model. The marcs model fails 
to reproduce the observations, regardless of the 
Sh used. 

4. Spatially-resolved line formation 

In Fig. |2]we show the distribution of the equiv- 
alent widths in the disk-centre granulation. The 
results from the 3D model have been con- 
volved with a point-spread function to account 
for the atmospheric smearing and the tele- 
scope's resolution, and with a Gaussian to ac- 
count for the spectrograph's instrumental pro- 
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Fig. 1. Centre-to-limb variation of the oxy- 
gen lines. Top: Oi 615.81 nm line (including 
blends). Middle: [Oi] 630.03 nm line (includ- 
ing blends). Bottom: Oi 777.41 nm line, only 
for the 3D model and for different scalings of 
the Hi collisions. 



file. The oxygen abundance was adjusted so 
that the predicted mean temporal and spatial 
equivalent widths matched the observed. 



For the O i 777 nm lines the NLTE results 
are shown, with S h = 1 . It is comforting to find 
that Sh = 1 gives the best agreement also here. 

5. Conclusion 

We looked at the centre-to-limb variation and 
spatially-resolved line formation of oxygen 
lines, comparing new observations with the 3D 
model. These lines are an important test of the 
model, because oxygen is very relevant in the 
recent revisions of the solar chemical compo- 
sition. 

Overall there is a very good agreement be- 
tween the predictions from the 3D model and 
the observations. Both in the centre-to-limb 
variation and the spatially-resolved line forma- 
tion for oxygen. These results give us confi- 
dence that the 3D model is realistic and appro- 
priate to derive the solar chemical composition, 
in particular the oxygen abundance. 

Acknowledgements. TMDP acknowledges financial 
support from Fundacao para a Ciencia e Tecnologia 
(reference number SFRH/BD/2 1888/2005) and 
from the USO-SP International Graduate School 
for Solar Physics under a Marie Curie Early 
Stage Training Fellowship (project MEST-CT-2005- 
020395) from the European Commission. This re- 
search has been partly funded by a grant from 
the Australian Research Council (DP0558836). The 
Swedish 1-m Solar Telescope is operated on the is- 
land of La Palma by the Institute for Solar Physics 
of the Royal Swedish Academy of Sciences in the 
Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofsica de Canarias. 

References 

Allende Prieto, C, Asplund, M., & Fabiani 
Bendicho, P. 2004, A&A, 423, 1 109 

Asplund, M. 2005, ARA&A, 43, 481 

Asplund, M., Carlsson, M., & Botnen, A. V. 
2003, A&A, 399, L31 

Asplund, M., Grevesse, N., & Sauval, A. J. 
2005, in Astronomical Society of the 
Pacific Conference Series, Vol. 336, Cosmic 
Abundances as Records of Stellar Evolution 
and Nucleosynthesis, ed. T. G. Barnes, III & 
F. N. Bash, 25-+ 

Asplund, M., Grevesse, N., Sauval, J., & Scott, 
P. 2009, ARA&A, 47, 481 



Pereira et al.: Testing 3D solar models against observations 285 



O I 777.19 nm O I 777.41 nm 1 777.53 nm 




| 0.3 r 

0.8 0.9 1.0 1.1 1.2 0.8 0.9 1.0 1.1 1.2 



Continuum intensity Continuum intensity 

Fig. 2. Distribution of equivalent widths over the solar granulation at disk-centre, for five oxygen 
lines. The 3D model results have been degraded to account for the atmospheric turbulence and 
the telescope's resolution. For each spatially-resolved spectrum the equivalent width has been 
computed and the resulting histogram of points is represented by the contours, as a function of 
the normalized local continuum intensity. The line profiles of the 1 777 nm lines were computed 
in NLTE, for Sn = 1. For the other two lines LTE was assumed. 
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